Deciding Unitary Equivalence Between Matrix Polynomials and Sets of Bipartite 

Quantum States 
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In this brief report, we consider the equivalence between two sets of m + 1 bipartite quantum states 
under local unitary transformations. For pure states, this problem corresponds to the matrix algebra 
question of whether two degree m matrix polynomials are unitarily equivalent; i.e. UAiV^ = Bi for 
< i < m where U and V are unitary and {Ai,Bi) are arbitrary pairs of rectangular matrices. We 
present a randomized polynomial-time algorithm that solves this problem with an arbitrarily high 
success probability and outputs transforming matrices U and V. 



INTRODUCTION 

With entanglement being one key component in the 
design and operation of quantum computers, it has be- 
come natural to treat entanglement as a resource which 
we extract from quantum systems and put to use. Under 
this interpretation, much research has been devoted to 
quantifying the amount of entanglement present in the 
state of a given system However, it was soon real- 
ized that no single quantification or entanglement mea- 
sure can fully capture a state's non-classical properties, 
and thus one must first stipulate a relative measure when 
asking how much entanglement some state possesses 
A common property of all meaningful measures is that 
entanglement between two subsystems cannot increase 
on average when manipulations are local, or applied to 
each subsystem distinctly; global actions are required to 
increase entanglement [3j. Because of the reversibility 
in unitary evolution, an immediate consequence of this 
is that for all entanglement measures, entanglement re- 
mains constant under local unitary operations (LU). As 
a result, studying LU equivalence is important since it 
identifies states that have the same amount of entangle- 
ment. 

With this motivation, we investigate the question of 
when two sets of bipartite states are simultaneously re- 
lated by a local unitary operation. More precisely, given 
two sets of states {po, p m } and {co, cr m } shared be- 
tween parties Alice and Bob, when is it possible for the 
duo to apply a fixed local unitary operation that pair- 
wise transforms pi ^> <7j for < i < to? Physically, this 
question arises when Alice and Bob are secretly given 
one of to + 1 possible initial states, and they wish to 
know whether, with certainty, it is possible to obtain one 
of to + 1 particular target states through local unitary 
operations alone. 

In the specific case of just a single pure state pair p and 
a, LU equivalence is decided by an equivalence in singu- 



lar values. Recently, significant progress has been made 
in the study of single-copy equivalence between multipar- 
tite pure states 0,0- Something not provided in these 
works, however, is a rigorous account of computational 
costs required to implement the described algorithms. In 
the case of generic bipartite mixed states, equivalence be- 
tween p and a is determined by a set of trace invariants 
@, 0], while the full solution to bipartite mixed state 
LU equivalence still remains open. The generalization 
of these questions to simultaneous LU equivalence be- 
tween multiple pairs of states has yet to be addressed, 
and such an investigation nicely complements previous 
work on simultaneous state transformations under global 
operations [8H10I] and simultaneous stochastic local state 
transformations between two pairs of pure states [llj . 

On its own, the above question generalizes the purely 
linear algebraic problem of deciding for to + 1 pairs of 
di x di matrices (Xj , Yi ) whether there exists unitary 
matrices U and V such that UXiV^ = Yi for all i. To 
our knowledge, this problem has not yet been studied 
either in the linear algebra community, although Radjavi 
has solved the special case of square matrices and U = V 
12j |. The problem can be phrased in a manner better 



suited for deeper analysis by introducing degree to matrix 
polynomials V{\) = £™ A**j and Q(A) = £™ A^. 
Two matrix polynomials are called unitarily equivalent 
if UVV^ = Q, and we see that UX.V^ = Y, for all i if 
and only if their corresponding matrix polynomials are 
unitarily equivalent. We also note a more general notion 
of matrix polynomial equivalence in which V ~ Q if there 
exists invertible constant matrices A and B such that 
AVB- 1 = Q. 

In this report, we present a randomized polynomial- 
time algorithm that decides whether two sets of bipar- 
tite pure states can be made equivalent by a fixed local 
unitary operation. For sets of iV-partite mixed states, 
the algorithm can be used to decide whether each pair 
is simultaneous equivalent under the same unilocal uni- 
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tary operation. These are special operations in which 
just a single party applies a local unitary while the other 
subsystems are left unperturbed. Our algorithm applies 
to sets of any size and the probability of failure can be 
made arbitrarily small since the randomness arises from a 
polynomial identity testing subroutine in the algorithm. 
The underlying technique of the algorithm also works to 
decide when two degree m matrix polynomials are equiv- 
alent in the more general sense of invertible transforming 
matrices A and B. 

Finally, we note that our result will also decide general 
(not just unilocal) LU equivalence of generic bipartite 
mixed states, although the set of "generic states" in our 
case is different than those in [1, 0] • Generic here means 
that the set of states to which our algorithm does not 
apply has measure zero. Specifically, our algorithm can 
be implemented on states that have distinct eigenvalues. 
If 9 = Ej Ci\fa)(fa\ with a > c l+1 and a = J2i C M)(<^I 
with > c' i+11 then p and a are LU iff c; = and 
\fa),\fai) are LU equivalent for all i. 

The actual problem we will consider is a bit more gen- 
eral than the one described in the previous paragraphs 
and will be called the Unitary Equivalence Problem 
(UEP): 

Suppose Gi and G 2 are sub-C-algebras of the 
rings C dlXdl and C d2Xd2 , respectively. For two 
sets of matrices {Ai}i = o,...,m and {!i}»=o,...,m with 
Xi,Yi £ C dlXd2 , decide if there exists a unitary 
solution U and V to the system of equations 



X 



= {UX i V* = Y i :U EGi,V eG 2 }. 



(1) 



The UEP formulation generalizes many different uni- 
tary equivalence problems. For instance, if we let G\ = 
C dlXdl and G 2 = C d2Xd2 , we recover the question of 
whether there exists general unitaries U and V such that 
UXiV' = Yi for all pairs (Xi, Yi). If we furthermore con- 
sider d\ = d 2 with one pair of matrices both being the 
identity matrix (7^ , 1^ ) , the question becomes whether 
UXiW = Yi for all i. An example of a nontrivial algebra 
Gi is the set {M ® I b : M € C axa } where ab = di. 

It is easy to see the connection between UEP and the 
simultaneous LU equivalence between bipartite states. 
The states of a di x ^-dimensional bipartite system 
can be represented as vectors \fa) in the product space 
C dl ® C da , and linear operators on this space correspond 
to physical actions on the system. By choosing some 
basis and \i) 2 for spaces C dl and C d2 respectively, 
any state can be written as \ip) — (I <g> where 
|$) = J2i=i N)iN)2- This allows for a bipartite pure 
state to be identified with the matrix ip 6 £dixd 2 so 
that the transformation — > (A®B)\tp) corresponds to 
ip — > AipB T . Consequently, simultaneous LU equivalence 
between states {\ipi)}i=o,„ m and {\fa)}i =o...m amounts to 
whether UipiV* 1 = fa for all i. For bipartite mixed states, 



the UEP is encountered only in the restricted setting of 
unilocal equivalence. Since mixed states themselves are 
represented by elements in C d i d 2xdid 2) un ilocal unitary 
equivalence between states p and a is the question of 
whether (U (g) Id 2 )p(U 1< ® Jd 3 ) = cr, which as noted above 
is an UEP instance. Note that in the case of simultane- 
ous unilocal equivalence of mixed states, the reduction 
to UEP applies to systems with an arbitrary number of 
parties. 



THE ALGORITHM 

As we will see in greater detail, the UEP can be solved 
by determining whether or not a particular system of 
quadratic equations has a nontrivial solution. One strat- 
egy sometimes helpful for dealing with quadratic con- 
straints is to relax the problem into a system of linear 
equations such that a solution to the new equations will 
solve the original with high probability. We demonstrate 
this idea on the problem of deciding whether two d\ x d 2 
(assume d 2 > d\) matrix polynomials V — J2"Lo 
and Q = YllLo are generally equivalent, i.e. V ~ Q. 
In other words, does the system of equations 



Xi = {AXiB- 1 =Yi : A € C dlXdl , 
B e C d2Xd2 ,0<i<m} 



(2) 



have a nonzero solution for invertible A and Bl Clearly 
Xi has such a solution iff there are nonzero invertible 
solutions to 



{AXi 
B € ' 



= Y t B : A e C dlXdl , 
Z d2Xd2 ,0 < i < m}. 



(3) 



There are (^(mrfj) linear equations in x'x which can be 
solved thus placing constraints on the 0(d 2 ) free vari- 
ables of A and B. A matrix solution space to x'i is then 
generated by expressing A © B in terms of the remain- 
ing free variables, and Xi has a solution iff there exists a 
nonsingular element in this space. 

A standard randomized algorithm for deciding whether 
a matrix subspace has a full rank element consists of eval- 
uating the degree O^) real polynomial \Det(A © B)\ 2 
for randomly selected values of the free variables. The 
Schwartz-Zippel Lemma states that for some n-variate 
polynomial f(xi, ■ ■ ■ ,x n ) over a field K and having de- 
gree no greater than d, if / is not identically zero, then 
Prob[/(a;' 1 , • • • , x' n ) = 0] < ryr where each x[ is indepen- 
dently sampled from some finite set Xcl 13- 1^. 

To use the Schwartz-Zippel Lemma for testing whether 
\Det(A © B)\ 2 is identically zero with success probabil- 

1d 2 

ity at least 1 — r^, one evaluates it on values randomly 
chosen from set IcR and decides a zero identity if and 
only if the evaluation output is zero. As any polynomial 
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number of linear equations can be solved in a polyno- 
mial amount of time in order to obtain the space A © B, 
we thus have an efficient method for deciding whether 
V ~ Q up to any probabilistic degree of certainty. We 
note that the Schwartz-Zippel technique can also be used 
in the study of bipartite entanglement distillation from a 
multipartite-party state 



16] 



To solve Xj we work analogously to xi but with addi- 
tional constraints enforced. Consider the system 



X ' = {AXi = YiB, XiB^ = A^Yi :A,A^eG u 

7-> nt , r. ^ ■ 1 



B, B^ e G-2; ,0 < i < to} 



(4) 



Then we have 



Theorem 1. x has a solution iff x' has an invertible 
solution A and B. 

Proof. If such a solution for x' exists, then A' AXi = 
XiB^B and AA^Y l = Y t BB^ . But these equations im- 
ply p(A^A)X l = X t p{B^B) and p(AA*<)Y t = Y iP {BB^) 
where p is any polynomial function. Let Xi denote 
the distinct eigenvalues from the combined spectrums 
X(A^A) U \{B^B). Let X be the Vandermonde matrix 
of the Xi, and v the column matrix whose entries are 
. Then the entries of X _1 v provide the coeffi- 



cients of a polynomial pit) such that p(A^A) = V At A 

and p{B^B) = y/WW . Note also that p(A^A) £ 
Gi and p(B^B) £ Gi- Define unitary matrices U = 

AVAlA^ 1 e Gi and V = bVWb' 1 £ G 2 . Then 
UXi = AXtx/WB^ 1 = YiBVbW 1 = YV. □ 

Matrix bases for G\ and Gi will contain no more than 
d\ elements so that x' represents 0(md\) linear con- 
straints on 0(d|) free variables. Indeed, two additional 
variable matrices M,N can be introduced to x' giving the 
equations AX t = YiB, XiN = MY h 5* = N, A f = M, 
A,M £ Gi, and B,N £ G2- A solution matrix space 
A B is generated, and like before, a polynomial iden- 
tity test can be applied to decide with arbitrarily high 
probability whether this space contains a nonsingular el- 
ement. If a nonsingular element is found, use the A and 
B to form unitaries U and V as in Theorem (JT]). 



CONCLUSION 

In this article we have studied the general problem of 
determining when a set of matrix transformations can 
be simultaneously achieved by a left and right unitary 



action. Physically, this corresponds to performing mul- 
tiple transformations between bipartite pure states with 
the same local action so that the amount of entangle- 
ment remains unchanged. Our analysis also extends to 
the situation of simultaneous unilocal unitary transfor- 
mations on TV-partite mixed states. We have developed 
a polynomial-time randomized algorithm that decides the 
problem with high probability and also provides a uni- 
tary solution if it exists. Related open questions con- 
cern simultaneous general LU equivalence between sets 
of bipartite mixed states and of states having more than 
three parties. However, for this latter question, it ap- 
pears the matrix polynomial and randomized techniques 
used above apply only to the types of LU equivalence 
considered in this report. 
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